Quantum criticality of the Lipkin-Meshkov-Glick Model in terms of fidelity susceptibility 
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We study the critical properties of the Lipkin-Meshkov-Glick Model in terms of the fidelity susceptibility. By 
using the Holstein-Primakoff transformation, we obtain explicitly the critical exponent of the fidelity suscep- 
tibility around the second-order quantum phase transition point. Our results provide a rare analytical case for 
the fidelity susceptibility in describing the universality class in quantum critical behavior. The different critical 
exponents in two phases are non-trivial results, indicating the fidelity susceptibility is not always extensive. 
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I. INTRODUCTION 

The Lipkin-Meshkov-Glick (LMG) model [1] was intro- 
duced in nuclear physics. It describes a cluster of mutually 
interacting spins in a transverse magnetic field. In condensed 
matter physics, this model is associated with a system of in- 
finite coordination number. In earlier time, scaling behaviors 
of critical observables have been studied by mean field anal- 
ysis (2J], while recently the finite-size scaling of this model 
was studied by the 1 /N expansion in the Holstein-Primakoff 
single boson representation J3[]and by the continuous unitary 
transformations (CUT) U ETSl- Meanwhile, a rich struc- 
ture of four different regions is revealed in the parameter 
space through a careful scrutiny on the spectrum [7]. Besides, 
the quantum criticality has been investigated by studying its 
entanglement properties H H OH E2l • Both the first- 
and second-order quantum phase transitions (QPTs) II 1311 have 
been revealed, in the antiferromagnetic and the ferromagnetic 
cases respectively 

Regarding the QPT itself, the ground state of a system 
would undergo a significant structural change at certain crit- 
ical point. This primary observation suggests a new descrip- 
tion of QPTs in terms of fidelity JH QlM [H M, El, Hfl 
I2ul^l2^l24ll25lp6ll . a concept introduced in quantum infor- 
mation theory 1131 . Mathematically it is the overlap between 
two ground states in which their driving parameters deviate 
by a small amount. However, the fidelity depends computa- 
tionally on an arbitrarily small yet finite change of the driving 
parameter. For this, Zanardi et. al. introduced the Riemannian 
metric tensor 11811. while You et. al. suggested the fidelity sus- 
ceptibility 11911 . both focus on the leading term of the fidelity, 
in order to explain singularities in QPTs. In addition, scal- 
ing analysis of these quantities has been informative: it helps 
understanding their divergence and the criticality of the sys- 
tem If2lll . and it also reveals the intrinsic relation between the 
critical exponent of some physical quantities and that of the 
fidelity susceptibility [22]. 

In this paper, we explicitly compute the ground-state fi- 



delity susceptibility and its critical exponent of the LMG 
model. Numerical analysis is also performed to check with 
our analytic calculations. We show that, the 1/7V expansion 
in the Holstein-Primakoff transformation is sufficient to de- 
termine the critical exponent of the fidelity susceptibility;^. 
In addition, we revealed two distinct critical exponents in two 
phases which is not a general feature. Therefore, our findings 
not only suggest another route on understanding the quantum 
criticality of the LMG model, but also show the fidelity sus- 
ceptibility is not always extensive in describing the universal- 
ity class of a quantum many-body system. 

This paper consists of five sections. In Sec. HJ we review 
the Hamiltonian, symmetry, and conserved quantities of the 
LMG model. In Sec. [Till we diagonalize the model Hamilto- 
nian and compute the fidelity susceptibility in the anisotropic 
model. In Sec. IIV1 we perform finite size scaling analysis 
and discuss the scaling relation between different exponents. 
Finally, we give a brief summary in Sec. [V] 



II. THE MODEL HAMILTONIAN 



The Hamiltonian of the LMG model reads 



~(s? + r S, 2 )-2W« + |(l + y). 
-■^(l+y)(S 2 -S 2 -A72)-2/^ 



(i) 

(2) 
(3) 



where <t k (k = x,y,z) are the Pauli matrices, S K = z^i<r l K l2, 
and S± = S x ± iS v . The prefactor 1/7V is necessary to en- 
sure finite energy per spin in the thermodynamic limit. It is 
understood that the total spin and the parity are the conserved 
quantities, i.e., 
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= 0. 



(4) 
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In addition, in the isotropic case y = 1, one has [H,S Z ] = 
and simultaneous eigenstates can be found. In the main con- 
text, the following parameter space is considered: A = 1, \y\ < 
1 , h > 0. We take h > as the spectrum is invariant under the 
transformation h «-» -/z. In addition, as a common practice 
we only consider the maximum spin sector S — N/2 which 
contains the lowest energy state. 



III. CRITICAL BEHAVIOR OF THE FIDELITY 
SUSCEPTIBILITY 

We briefly review of the concept of the fidelity susceptibil- 
ity here. Suppose there is a Hamiltonian of a general form 



H — H()(y) + f(h)H; 



(5) 



for Hi is defined as the driving term of the system, which sim- 
ply does not commute with Hq. The function /(/;) coupled to 
Hj is often considered as the linear external fi eld f (h) - h. 
Then the fidelity susceptibility is defined as lfl8i Il9ll 



df(h) 



dh 



(6) 



where E„ and \n) stand for the « th eigenenergies and eigen- 
states of the (whole) Hamiltonian respectively. 

The fidelity susceptibility is well-defined for a non- 
degenerate ground state of the continuous variable h, but it is 
not suitable to deal with states with good quantum numbers. 
The LMG model undergoes ground state level crossing when 
y = 1, the ground states are assigned the magnetization as the 
quantum numbers. 

We put our focus on the fidelity susceptibility for an arbi- 
trary isotropy M < 1. One resolution is to use the Bethe- 
Ansatz solution [|27il28ll . which is rather complicated. So we 
adopt the 1 /N expansion method which was used extensively 
by Dusuel and Vidal Jj, |5|], that corresponds to the large N 
limit. 

The l/N expansion method is done under the Holstein- 
Primakoff boson representation [3] framework. In low en- 
ergy spectrum the spin operators in the S = N/2 subspace 
are mapped into boson operators: 

S z — S - a'a, 

S + = (2S - a'a) ll2 a = N l/2 (1 - cJa/N) ll2 a = S'l, (7) 

where a{a ") is the standard bosonic annihilation (creation) op- 
erator satisfying [a, a 1 "] = 1. The above transformation is 
valid when h > 1, but when < h < 1 it can also be used 
through semi-classical treatment This representation is 

also known as the spin-wave theory. It is well adapted to the 
computation of the low-energy physics when (a^a)/N <K 1. 
After inserting these expressions of the spin operators in Eq. 
Q, one can approximate the square roots as one and express 
the result in normal ordered form with respect to the boson 
vacuum state. Keeping terms of order (1/N)~ l , (1/N)~ 1 ^ 2 and 



(l/N)° for h > 1 (in which the approximation is justified), the 
Hamiltonian becomes 



H = -hN + (2h - 1 + y)a'a - (a' 1 + a 2 ) . (8) 



The above Hamiltonian can be diagonalized by a standard Bo- 
goliubov transformation 

a 1 " = cosh(0/2)fo' + sinh(0/2)£, (9) 
a = sinht©^)/? 1 " + cosh(0/2)£, (10) 

where /?(£>' ) is the quasi-bosonic annihilation (creation) oper- 
ator, and 



tanh[0(/i > 1)] 



l-y 



2h-\+y 
then the Hamiltonian is diagonalized as 



(11) 



H = -h(N + 1) + 2 y/(h - \)(h - y) [b^b + X - j . (12) 

Thus the low-energy spectrum of the model is mapped to the 
spectrum of a simple harmonic oscillator. The eigenstates are 
just {|n)}, where b^b\n) = n\n). We consider the driving Hamil- 
tonian Hi responsible for the QPT, 



-2S,. 



(13) 



By transforming them into combinations of b and £>' opera- 
tors, the fidelity susceptibility is calculated as 



d-r) 2 



32(A- lf(h-y) 2 



(14) 



The derivation above is only valid for h > 1, for < h < 1 
the calculation is actually similar to the above case of h > 1, 
provided that one first rotates the z axis to bring it along the 
classical spin direction. We do not show it explicitly here, but 
interested readers are recommended to refer to Ref. |4j, |5[]. 
We simply quote the main result, after all the procedures the 
Hamiltonian becomes: 



+ 2 ^(1 -#)(!- y) (fctft + ij(.15) 



(1+h 2 ) l-y 
H = -N - 

2 2 

The driving Hamiltonians also takes a different form: 

i 

= -2(- Vl -h 2 S~ x + hS z ), (16) 

for the HP transformation is done on the S operators. The 
fidelity susceptibilities are then obtained accordingly: 



N 



X F = 



h 2 (h 2 - yf 
+ ^ — (17) 

. . . 22 - V LI) 



4V(l-A 2 )(l-y) 32(l- r ) 2 (l-/* 2 ) 
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FIG. 1 : (color online) The fidelity susceptibility in response to h as 
a function of h at y = 0.5. The inset denotes the difference between 
the fidelity susceptibility and the extensive term in Eq. J17b . 



Thus we obtained \ F of the anisotropic LMG model in large 
N limit. We first see the effect of isotropy to the fidelity sus- 
ceptibility. It dominates when h < 1, but fades out for large 
h. Especially in the isotropic limit, when y — > l,x F diverges 
when h < 1, but tends to zero when h > 1. This is the ef- 
fect of the level-crossing points in the thermodynamic limit. 
They together form a region of criticality, and the system un- 
dergoes continuous level crossing. The fidelity susceptibility 
responds drastically while moving along h. But when h > 1, 
there are no further critical points, x F naturally measures zero 
when moving along h because we have [Ho, Hj] = 0. 

An interesting observation is x F behaves extensively when 
h < 1 even in the large N limit. When discarding the exten- 
sive part of Eq. ( fTTT i. we arrive a zero point at h = y/y, which 
does not fit with numerical analysis [Fig. Q~). This discrep- 
ancy may be eliminated by adopting other transformations of 
the driving Hamiltonian. Particularly, the flow of operators in 
the LMG model haven been studied by the continuous uni- 
tary transformation (CUT) method [4, 5fl. However, such dis- 
crepancy would not hinder us from getting the correct critical 
exponent of the fidelity susceptibility. 

Let us emphasize the intensive property of the fidelity sus- 
ceptibility, which measures the average response to some driv- 
ing Hamiltonians. Its divergence should correspond to a criti- 
cal point of a second-order QPT rather than to the increasing 
system size. In order to predict the critical exponent correctly, 
we should average the fidelity susceptibility whenever neces- 
sary. To the leading order, Eq. ( fT7b becomes 



N 



1 



4V(l-/* 2 )(l-y) 



(18) 



Then it comes to a key result of our paper: x F bears differ- 
ent critical exponents across the critical point. It diverges as 

i ~ 
(1 - h)~- when «<!,(«— 1) when h > 1. It is unlike the 

Ising model in a transverse field lfl5ll nor the one-dimensional 
asymmetric Hubbard model Ir22tl . where the critical exponent 
is a single number over the phases. 



IV. FINITE SIZE SCALING ANALYSIS 

To illustrate the scaling behavior of the fidelity susceptibil- 
ity, we perform the exact diagonalization (ED) to solve the 
spectrum of H and then calculate the corresponding fidelity 
susceptibility numerically. 

Let us recall the fidelity susceptibility scaling analysis per- 
formed in the asymmetric Hubbard model M22I1 . According to 
the scaling ansatz l29ll and the obvious power-law divergence 
observed in Fig. Q] the rescaled fidelity susceptibility around 
its maximum point at /i max is a simple function of a scaling 
variable, i.e. 



X, 



■X F 



x,. 



= f[N v (h-h max )], 



(19) 



where f(x) is the scaling function and v is the correlation 
length critical exponent. This function is universal and does 
not depend on the system size, as shown in Fig. [2]for cases of 
y = 0.5,0 and y = -0.5. Remarkably, the critical exponent v 
for three cases are very close. This observation strongly im- 
plies that v is a universal constant and does not depend on the 
parameters y and h. 

In recent studies of the fidelity susceptibility in critical phe- 
nomena, it was pointed out that the intensive fidelity suscepti- 
bility scales generally like If2lll22tl 



1 



X, 



\h-h c \ a 



(20) 



around the critical point. In the last section, we have already 
obtained 



a = 



h > 1 
< h < 1 



(21) 



which is also a universal constant. Then if the maximum point 
of the intensive fidelity susceptibility scales like 



(22) 



the scaling ansatz also implies another important relation, i.e. 



A* 

a — —. 

v 



(23) 



We try to confirm this equality in numerically. In Fig. f2] 
Eq. ( fT9] > is best fitted with v ^ 0.665. The case to determine 
// is more subtle. It is because Eq. ( fT9] > remains the same 
form even for averaged x F , but the maximum of^ f does not. 
To resolve this problem, we first determine [i from the "bare" 
X F . By using least square fit method, we evaluated "bare" // 
for at different y. The numerical details are shown in table 
U However, the exponent fi does not converge perfectly. We 
compare the /u obtained in a range of [2 12 , 2 ], and those from 
the range [2 8 , 2 16 ]. The results converge better for larger scal- 
ing regions. According to the trend of fj. in larger system sizes, 
we roughly estimate fx = 1.33 with three effective digits [Fig. 

When h > l,x P is observed to be intensive [Fig. [T|. With 
the estimated fi and v, the equality d23l is satisfied with a — 2. 




FIG. 2: (color online) The finite size scaling analysis is performed for the case of power-law divergence at y = 0.5 (LEFT), y = (MIDDLE) 
and 7 = -0.5 (RIGHT) for system sizes N = 2"(n = 12, 13, 14, 15, 16). The fidelity susceptibility is considered as a function of system size 
and driving parameter is a function of N v (h - /i max ) only, with the correlation length critical exponent v = 0.665. 



TABLE I: Scaling exponent ji at various y, obtained by sampling system size in different range. 



r 


0. 


i 


0.5 


0.2 







-0.2 


-0.5 


H(N e [2 s , 2 16 ]) 


1.3221 ± 


0.0006 


1.3264 ±0.0004 


1.3267 ±0.0004 


1.3280 


± 0.0004 


1.3283 ± 0.0003 


1.3285 ±0.0003 


H(Ne [2 12 ,2 16 ]) 


1.3250 ± 


0.0003 


1.3285 ±0.0004 


1.3295 ±0.0002 


1.3299 


± 0.0002 


1.3302 ± 0.0001 


1.3304 ±0.0001 



M y=0.8 




qI 1 1 1 1 1 1 

10 12 14 16 

log 2 iV 



FIG. 3: (color online) The finite size scaling is performed for the 
maximum of the fidelity susceptibility. 



On the other hand, when h < l,y F IN is the intensive quantity. 
For^ p oc N M , 

oc A^-D. (24) 

N 

Thus n =i 0.33, this will give the relation a — 1/2. These two 
values of a are consistent with our analytic calculation in the 
last section. 

The exponent fi, v can also be discussed from the scaling 
ansatz at the critical point rather than the maximum point 
of a finite system, as shown by Vidal, Dusuel, and Barthel 
Based on their approach, the critical exponent v takes 
the value of 1/3, and is independent of magnitude of y. Our 
results on the maximum simply agree with this value and can 
be generalized to other models where the precise critical point 



is not known. 

Another scaling analysis is to examine how /z max tends to 
the critical point h c — 1 . It should scale like 

h c -h mllx KN- d , (25) 

in the large N limit. We find 6 - 0.66 with two effective digits 
for various y. 

In short, we can confirm that the exponents fi, v, and 5 of 
the fidelity susceptibility do not depend on the value of y and 
h. They are universal constants for the LMG model and are 
related to the critical exponent of the fidelity susceptibility a. 

V. SUMMARY 

In summary, we computed explicitly the fidelity suscepti- 
bility and its critical exponent of the LMG model at different 
isotropy. We confirmed the different critical exponents in two 
phases numerically by ED, which is a rather non-trivial result. 
Several scaling exponents are also found in consistence with 
previous studies. Since the fidelity susceptibility is believed to 
be able to characterize the universality class of quantum phe- 
nomena, our results therefore provide a rare explicit case for 
the study of fidelity susceptibility. 
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